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Qh! This is the continuation of a previous article, in which the Bjorken and 



Abstract 



Voloshin sum rules were interpreted as statements of conservation of prob- 
ability and energy. Here the formalism is extended to higher moments of 



X 

5^ , the Hamiltonian operator. From the conservation of the second moment of 

the Hamiltonian operator one can derive a sum rule which, in the small veloc- 
ity limit, reduces to the Bigi-Grozin-Shifman-Uraltsev-Vainshtein (BGSUV) 
sum rule. On the other hand, the conservation of the third moment of the 
Hamiltonian operator gives a new sum rule, which is related to the matrix 
element of the heavy quark counterpart of the Darwin term in atomic physics. 
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The general case of the higher order moments is also discussed. 



Typeset using REVT^ 



I. INTRODUCTION 



In a previous article |]T|, the relationship between various conservation laws and various 
sum rules has been explored. In the heavy quark limit, the hadronic resonances appear as 
eigenstates of the light degrees of freedom under the static color field of the heavy quark. 
Under a 6 — c decay, the Isgur-Wise form factors fn'^w) are defined to be, up to a kinematic 
factor, the overlap of the initial ground state brown muck |0) moving with velocity v and 
the final (ground state or excited) brown muck \n') with velocity v', 

cpniw) = (n'lO), (1.1) 

with w = V ■ v'. By the completeness of the set of eigenstates, we have the master sum rule 

E(0|XK)(n'|0) = (0|X|0), (1.2) 

n' 

where X is an arbitrary operator. The case of interest is when X is conserved, i.e., H' and 
X commute and the |n')'s are eigenstates of X. Then the left hand side of Eq. ( |1.2|) is just 
the sum of the squares of the Isgur-Wise form factors f{w) weighted by their respective 
eigenvalues of X. For X = U, the identity matrix, this master sum rule is equivalent to the 
well-known Bjorken sum rule pj-Q, 

Y^\^n'{w)\' = l. (1.3) 



can 



On the other hand, with X = H', the hamiltonian operator in the v' frame, Eq. ( |1.2|) 
be reduced to the Voloshin sum rule (which will be discussed below). Lastly, a new sum 
rule can be obtained by choosing X = P', the parity operator in the v' frame. This new 
"parity sum rule" can lead to model independent lower bounds of (po'{w), the ground state 
Isgur-Wise form factor. 

In this paper, we are going to consider the case X = H''' with k > 2. When X = H'^, 
the sum rule obtained is closely related to the Bigi-Grozin-Shifman-Uraltsev-Vainshtein 
(BGSUV) sum rule derived in Ref. When X = H'^, a new sum rule, relating the Isgur- 
Wise form factors and the matrix element of the heavy quark counterpart of the Darwin 



term in atomic physics, is obtained. The general case of X = H with > 4 is briefly 
discussed. 

II. SUM RULES ON H''^ 

By considering the operator product expansion of heavy quark currents in the context of 
QCD sum rules, the authors of Ref. have been able to obtain several sum rules relating 
the Isgur-Wise form factors in the small velocity limit. Besides the Bjorken and Voloshin 
sum rules, they obtained the new BGSUV sum rule, 

T.Mw)fE'^, = lf^liv-v'y, (2.1) 

n' 

where En' = m^n' — ttib is the mass difference over the ground state meson, and 

= {BWDfh\B). (2.2) 

(Notice that the present notations are different from those in Ref. 0. In particular the 
normalizations of the meson states differ by a factor of \/2Mb-) 

Physically /i^ denotes the square of the heavy quark momentum inside the B meson, 
which is also the square of the brown muck momentum in the hadron rest frame. 

^^l = {0\p'\0). (2.3) 

Evidently this sum rule expresses conservation of the second moment of the energy operator, 
and we will see that this sum rule indeed can be reproduced as a special case of our master 
sum rule. 

Putting X = H'^ in the master sum rule, we have 

j:\^Ay^)\'Aml = m"\0), (2.4) 

n' 

with Arrin' = m-y-n' —rrib as the mass of the "brown muck". To see how the quantity {0\H'^\0) 

b 

can be evaluated, it is worthwhile to review briefly the derivation of the Voloshin sum rule 
The Voloshin sum rule is obtained by putting X = in Eq. (|1.2| ). 
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^Am„,|y;,,H|2 = (0|i/'|0). (2.5) 

n' 

To calculate the right-hand side, note that |0) and |0') are related by a Lorentz transforma- 
tion L which boosts from the v' frame to the v frame. 

|0) = L|0'). (2.6) 

Then 

(0|i/'|0) = (0'|L-^i7'L|0'). (2.7) 

The Lorentz transformation dictates that the hamiltonian is transformed as the time com- 
ponent of a four vector. 

L-^H'L = w{H' + {v-v')-f). (2.8) 

Since the ground state |0') is at rest in the v' frame, {0'\p'\0') = 0. On the other hand, 
{0'\H'\0') = A = niB — mb is the mass of the brown muck. Hence, 

{0\H'\0) = w{0'\H'\0') + w{v-v') ■ (O'Ip'IO') = wA, (2.9) 

as conjectured in Ref. [^. Together with Eq. (P75|), the Voloshin sum rule is reproduced: 

X]Am„/ !¥?„/(«;) 1^ = wA. (2.10) 

n' 

The argument above can be generalized to evaluate (0|if'^|0). 

{0\H'^\0) = w^{0'\H'^\0') + w\v - v'Uv - v')j{0'\p',p'^\0'). (2.11) 
with the terms proportional to {0'\p'\0') vanishing. By the relation 



{0'\p[p',\0') = mp"\0')S^, = (2.12) 



we end up with 



j:\^Aw)\'Aml, = {0\H"\0) 

n' 

= w^A^ + \w^fil{v-v'f. (2.13) 



(Strictly speaking, the matrix element ( |2.12|) contains in general also an antisymmetric part 



in the indices However, for the purpose of deriving the sum rule ( |2.13| ) this part can be 



neglected, as it vanishes under multiplication with {v — v')i{v — v')j in Eq. ( |2.11| ).) 

Using Eqs. (pTBI), ( |2.10| ), ( p.l3[ ), and the relation En' = Am„/ — A, the left-hand side of 



Eq. ( |2.1| ) can be evaluated. 

^ \^n'{w)\^El, = \^niw)\\A^ - 2AAm„, + Ami) 

n' n' 

= A\w-iy + lw^ljl{v-v'f. (2.14) 
So in the small velocity limit, i.e., ~ 1, 

J:Mw)\'EI, = 1^,1{v-v')\ (2.15) 

n' 

and the BGSUV sum rule is reproduced. 

Noting that {v — v'Y = — 1, Eq. ( p.l4|) can be formally expanded about w = 1 as 



Y.\^A^)\''El = l^,l{w-l) + .... (2.16) 

n' 

On the other hand, we know that in general, for a state with orbital angular momentum 
/ > 0, \ipn'{w)\'^ ^ {w — 1)'. For / = 0, \ipn'iw)\^ ~ (w — 1)° for the ground state and 
~ (w — 1)^ for excited states. This determines the behavior of the yj's near the point of zero 
recoil. 

The zeroth order expansion of Eq. ( |2.16| ) about the point of zero recoil is the trivial 
identity = 0, while the first non-trivial term is 

5: \^A^)\'El, = |^,l{w - 1) = If^liv - v')\ (2.17) 
1=1 

where v and v' are the four velocities of the initial and final heavy quarks. Then the above 
modification of the BGSUV sum rule is saturated by the P-wave states. (In contrast, the 
original BGSUV sum rule runs over all excited states.) 

In the meson sector, this modified BGSUV sum rule becomes 



E2 {Effff \r'^,{l)? + E4 1^1(1)1' = (2-18) 

q r 

where the P-wave Isgur-Wise form factors r(w)'s are defined in Ref. Although we have 
derived this modified BGSUV sum rule in the context of B decays, the same discussion holds 
for Ab decays as well, for which the sum rule reads 

T.^ {E['^y \a^^\l)\' = Ifil, (2.19) 

where /i^ is defined in analogy to /i^ 

fil = {Ab\^^Dfb\At,), (2.20) 

and the P-wave Isgur-Wise form factors cr(w)'s are defined in Ref. 0. 



III. SUM RULE ON H'^ 

In this section we will derive a sum rule from the conservation of the third moment of the 
hamiltonian. Bearing in mind that H and p are operators which in general do not commute, 
putting X = H'^ in the master sum rule gives 

El^„,(^)PA<' = (0|^10) 

n' 

= w''{0'\iH' + {v-v') ■ p'){H' + {v- v') ■ p'){H' + {v- v') ■ p')\Q') 
= w^{0'\H'^\0')+w%v- v'),{v- v')j{0'\{H'p',p'^ + p[H'p'. + p',p'^H')\0') 
= w'^A^ + w^AfiKv - v'f + \w^iil{v - v'f, (3.1) 

where 

/ij = {Q\{piHpi-piPiH)\Q) = (0|p,[i/,pi]|0). (3.2) 
E \^n'{w)\^El, = E \y^niw)W-A^ + SA^Anin' - 3AAm2, + Ami,) 

n' n' 

= -A^ + SA'^w - 3A{w^A^ + Iw^filiv - t;')') 

+ {w^A^ + w'^Afiliv - v'f + Iw^'filiv - v')^) 
= A^w- If + w\w - l)Afil{v - try + Iw^fi^iv - v')\ (3.3) 



Then 



The first, second and third terms are of order {w — 1)^, {w — 1)^ and {w — 1)^ respectively. 
As a result, when expanded about the point of zero recoil, we get 

El^n'HP^^' = M(^-l)- (3.4) 

1=1 

In the heavy mass limit, the commutators of the momentum and hamiltonian brown 
muck appearing in this matrix element can be written only in terms of heavy quark and 
gluon fields. The result for fi^ reads (details of the derivation can be found in the Appendix) 



f^l = -{B\Q^'giD'^F,^^rQ\B). (3.5) 



A sum rule similar to Eq. ( |3.4|) is obtained in the framework of [||, from considering the 
third momentum of the hadronic tensor (in the terminology of [P], the "fourth sum rule"). 
We have explicitly checked that the two methods lead to identical results. 



The matrix element /x?^ (|3.2|) has a special significance for heavy hadron physics. The 



dimension-6 operator in (|3.5| ) appears in the heavy quark effective theory lagrangian at 
order In nonrelativistic language (in the rest frame of the heavy hadron), fi^ is the 

expectation value of the Darwin term in the Pauli equation satisfied by the heavy quark field 
Q. A similar quantity has been defined as pfj in Eq. (27) of ^ where its contribution to the 
mass of a heavy hadron has been computed (see also llO|). Also, /i^ is directly connected to 
the third moment of the distribution function F{x) for inclusive semileptonic decays of heavy 
hadrons into final massless quarks. Its value determines the asymmetry of F{x) between 
positive and negative values of the scaling variable x. Finally, a quantity related to /i^ is 
needed when computing mass corrections of order to the total semileptonic width of 

a heavy hadron ||^ . 

Usually the matrix element ( p.5| ) is evaluated with the help of the factorization approxi- 
mation P,p!0|. For this the equation of motion of the gluon field is used D^Fq^ = g J2q Qlof^Q, 
where the sum runs over all light quarks. Then the Fierz identity is applied followed by the 
insertion of the vacuum state. The result reads 



where the case of a B meson has been considered. To turn this into a numerical prediction, 
a scale must be chosen at which the coupling is to be evaluated. We employ for this a low 
scale of the order = 0.5 GeV, where = 0.31 [H]. For we use an average value of the 



lattice QCD results fs = 175 ± 25 MeV [O. Together with the B meson mass tub = 5.28 



GeV these values give the following estimate 



/i^ = 0.140 ±0.043 GeV^ (3.7) 



On the other hand, the sum rule ( |3.4|) allows a direct evaluation of this quantity, in 
terms of the contributions of the excited states. Compared with the factorization approach, 
this method has the advantage of being free of uncertainties connected with the value of 
the factorization scale. Therefore it provides us with an independent test of the vacuum 
insertion approximation for hadrons containing one heavy quark. 

Expressed in terms of the contributions of the P-wave mesons, the sum rule ( p.4|) takes 
a similar form to ( p.l8[ ): 



2Ei/,\n/,{l)\' + 4El,\Ty,{l)\' = 1/^3 . (3.8) 

We have neglected here the contributions from the continuum and the higher radial excita- 
tions. The contributions on the left-hand side can be written [^], via the Bjorken sum rule, 
in terms of the slope p of the Isgur-Wise function at the equal- velocity point v = v': 

aw) = l-p\w-l) + --- . (3.9) 

One has 

/ = i+|n/2(l)P + 2|r3/2(l)P, (3.10) 

where again we have assumed saturation by the lowest-lying P-wave states only. Combining 
]8|) and (|3.10|) , the following bounds are obtained for /x^ 



^E!/, - i) < /ij < 3^3% (/ - i) . (3.11) 
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The excitation energy of the sY = | P-wave multiplet can be directly extracted from 
experiment. In the charm sector the masses the two members of the doublet are m£,^ = 
2423 ±3 MeV and rriD* = 2464 ±4 MeV (average values over the isospin doublet). These are 
expected to have arisen from a degenerated multiplet in the infinite mass limit, with a mass 
mpg^j — l/8(3?Tj,£i^ +5m£,*) = 2449 ±4 MeV. The average mass of the corresponding S-wave 
charmed states is mg = l/4(m/) + 3m/).). Using = 1867 MeV and m£, = 2008 MeV 
gives ms = 1973 MeV, from which the excitation energy E3/2 = 476 ± 4 MeV is obtained. 

Recently, the discovery of the sj'' = I"*" family of P-wave bottom mesons has been 
announced |p!3| , p!4| . The first of these collaborations quotes the following mass values: m^^ = 
5725 MeV and m^* = 5737 MeV, which combine to an average mp^^^ = 5732.5 MeV. The 
average mass of the S-wave bottom mesons is ms = 5279 MeV, which gives an excitation 
energy £'3/2 = 453.5 MeV. The difference from the charmed mesons' case can be attributed 
to spin-independent corrections proportional to l/m^fe. We will use in the following this 
latter number, which is presumably closer to the static value of this quantity in the infinite 
mass limit. 

None of the s^^ = I"*" P-wave states has been observed and the only information we have 



about them comes from model calculations [jl5|,|T6[. Their masses in the charm sector are 



expected to lie in the range mp^^^ = 2300 — 2400 MeV, from which an excitation energy 
of about Ei/2 = 380 ± 50 MeV can be extracted. In analogy with the preceding case one 
expects this value to decrease slightly when the mass of the heavy quark is increased. Thus 
we will use in the following E1/2 = 360 ± 50 MeV. 

The other parameter entering ( |3.11| ) is the slope of the Isgur-Wise function ai w = 1, 



which has been measured by the CLEO collaboration, with the result [lT7| , p!8[| = 0.84 ± 
0.12 ± 0.08. Combining these parameters, the bounds (|3.11| ) read 



0.083t°:°^^ GeV^ < 4 ^ 0.1651°;°^° GeV=^ . (3.12) 

Comparing with ( ^.7| ) one can see that the factorization approximation gives a rather 
j;ood estimate of the matrix element (|3.5|) . Despite the large errors on the lower limit in 
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( P.12| ), it seems justified to conclude that, at the very least, the factorization approximation 
provides us with the correct sign of this quantity^. 



IV. HIGHER MOMENTS OF H' 

Sum rules similar to the ones discussed in the preceding sections can be derived for even 
higher moments of the Hamiltonian operator. The general form for these sum rules can be 
obtained by using X = {H' — A)^ in the master sum rule (|1.2| ) with k > A. This gives 

5:|^„,H|2e^, = (0|(/7'-A)» = {0'\{wH' + w{v-v')-p'-Af\0'). (4.1) 

n' 

The matrix element on the right-hand side can be expanded in powers of v — v' and only 
the terms of even order give a nonvanishing contribution. Of these, only those quadratic in 
V — v' will contribute to order {w — 1) (since {v — v'Y = w'^ — 1). One obtains in this way 

E Mw)\'Et = w\v-v'Uv~ v%{0'\p[iwH' - A)'''p'^\0') + 0{{w - if) (4.2) 

n' 

= ^w^w' - lW\p[{wH' - A)^-V:|0') + 0{{w - If) 

= ^{w - i){o'\pkh' - A)'=-v:io') + oiiw - m . 

The matrix element on the right-hand side can be written as (removing the primes for 
clarity) 

ftl = {0\p,iH - A)'=-V.|0) = mp^, H]{H - A)'^-'[H, p,]|0) (4.3) 

fc-4 

= (0|[p„ H][Hy^., [H, p^]]■■■m■ 
For k = 2,3 the corresponding matrix elements have been denoted as /x^ and /i^, respectively. 



^For the purpose of comparison, we quote also the result of a calculation based on the ISGW model 
iQ] which gives ti/2(1) = ''"3/2(1) = 0.310. Together with the excitation energies given in the text, 
one obtains from (^.8]) = 0.067 GeV^. The slope of the Isgur-Wise function comes out equal to 



0.538, well under the CLEG result |17,18|, which seems to imply that /x^ is also underestimated. 
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The multiple commutator in ([4.3|) can be computed explicitly with the help of ([X9|). For 
example, 

[H, [H, p,]] = [Ht - H,, [H, p,]] = p,] - [H,, [H, p,]] , (4.4) 

where the notations are explained in the Appendix. The first term is equal to 

- tj^[H, p,] = Jdx (^Q\-t D^)igF^,Q + Q^igF^,{-iD^)Q + g{DA)Q] , (4.5) 



whereas the second one can be explicitly evaluated with the help of (^). By making use 



of the equation of motion (^) one finds that the latter just cancels the first two terms in 
LSI) so that one obtains 



[H, [H, pi\] = I dxQ^g{D,F,,)Q , (4.6) 

or, more generally, 

k-4: 

\lC~^[H, Pi] ■■■] = Hf-' I dxQ^g{D',-'F,,)Q . (4.7) 

Finally, one has for /i^ ( |4.3|) 

fit = -I J dxQ^gFoiQ ■ i-tf-' J dyQ^g{Dl-^F,i)Q (4.8) 
= {-if [ dxQ^g^FM-'Fo^)Q . 



In the last step we have made use of the fact that the quark Q is infinitely heavy and 
thus heavy quark pair creation is rigorously forbidden. This implies that a first-quantized 
(quantum-mechanical) description for the heavy quark is exact in this limit, with operators 
satisfying simple multiplication rules. The final result in ( |4.8| ) has been translated back in 
a second-quantized language^. 

The general form of our sum rule is thus 



similar argument has been used implicitly in the derivation of the BGSUV sum rule in section 
II when identifying the operator piPi in ( ^.3| ) with h[iD)'^h in (^.2|). 
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(4.9) 



1=1 



with only P-wave states contributing on the left-hand side. The matrix element appearing in 
this sum rule involves (for k> A) the chromoelectric field and its time derivatives. The same 
quantities have been shown in to play an important role in the inclusive decays of heavy 
hadrons to states containing a final massive quark. More precisely, the matrix element is 
connected to the (fc — 2)*'' moment of the so-called "temporal" distribution function for these 
decays. Therefore the most important application of the sum rules (|4.9| ) can be expected to 
be constraining the values of these matrix elements and thereby giving useful information 
about the temporal distribution function. Also, the operators whose matrix elements define 
/ifc might appear in physical applications as power corrections to the infinite quark mass 
limit. However, from dimensional arguments they are expected to appear suppressed by a 
factor of 1/Mq~^ and hence be unimportant when k is large. 



CK.C would like to thank Mark Wise for valuable discussions. His work was supported 
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This appendix contains the technical details of the computation of the matrix element 
of brown muck operators which appears in the definition of /i^ ( p.2|) . Our final aim is to 
express it only in terms of matrix elements of the heavy quark and gluon fields. We begin by 
defining the momentum and hamiltonian operators for the total system (brown muck plus 
heavy quark): 
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APPENDIX: 



Pt=P + Ph 



(Al) 



Hi 



t — 



(A2) 
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We will only need in the following the explicit form of the heavy quark operators: 



Ph 



J dxQ^{-iV - gA^t^)Q (A3) 
Hh = J dxQ^a ■ (-zV - gA'^e) + l3m]Q . (A4) 

Q is the heavy quark field operator and are the generators of the SUc(3) gauge group, 
satisfying the commutation relations [t", t^] = ifabcf^- The equation of motion satisfied by 
Q is 

iDoQ = [a ■ {-iV - gA'^e) + f3m]Q (A5) 

with Df, = d^ + igAp^ 

Let us consider the commutator between the heavy quark hamiltonian and the heavy 
quark momentum. It can be written as 

[Hh, ph^ = Jdx {Q^[a ■ i-tV - gA'^t'') + f3m]{-tW - gA)Q 
-Q\-iV - gA)[a-{-iV - gA''t'')+ (3m]Q] 
= jdx ^Q\-i Do){-tV - gA)Q - Q^{-tV - gA){tDo)Q^ (A6) 
= Jdx {Q^itDo, i-tV - gA)]Q - tdo[Q^i-tV - gA)Q]} 

= JdxQUgF,TQ-^^- 

In the second equality the equation of motion for the field Q ( |A5| ) has been used. The 
field tensor is defined as F^^ = d^A^ — duA'^ — gfabcA^^Al . The time- derivative on the 
right-hand side can be computed with the help of the equation of motion for Heisenberg 
operators idA/dt = [A, H], with H the total hamiltonian (|A2|) . In this way one obtains the 
commutator between the brown muck hamiltonian and the heavy quark momentum 

[H, = -J dxQhgFSrQ . (A7) 

So far everything has been completely general and we have not made use of the fact 
that the quark Q is heavy. The important simplification which appears in the heavy mass 



14 



limit is that it makes sense to speak about the eigenstates of the brown muck in a hadron 
moving with a well-defined momentum [^. This means that in this limit the brown muck 
hamiltonian commutes with the total momentum: 

[H,p,] = [H,p] + [H,ph]=0. (A8) 



From this relation and ( |A7D , one immediately obtains the needed commutator between the 
brown muck operators 

[H,p,] = J dxQUgF,TQ- (A9) 
It is now a simple matter to compute fi^ (|3.2|). It can be written as a double commutator 

/x^ = l(0|[p„ [H,pm. (AlO) 

where the relation 

{0\[H,p']\0) = {0\p,[H, p,]\0) + {0\[H, p,]p,\0) = (All) 

has been used. The vanishing of this matrix element is due to the fact that the states are 
eigenstates of H. Further, note that the first pi can be replaced in ( |A1U| ) by —phi since the 
external state (the heavy hadron is assumed at rest) is an eigenstate of the total momentum 

Pt\0) = 0. (A12) 

Thus one obtains 

/^x = -^(o|[p'^- [H, p^]m 

= -^(0| [/ d'xQ\-zd, + gAinQ, J d'yQHgF,TQm 
= i(0| Jd'xQ^gD^FS^eQlO) 

= ^(0| J d'xQ^''g{D^F,^^)t''Q\Q) , (A13) 
which is the result quoted in the text (|3.5|). 
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